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Abstract 

We consider the lagrangian of a self-interacting complex scalar field admitting generically 
Q-balls solutions. This model is extended by minimal coupling to electromagnetism and to 
■ gravity. A stationnary, axially-symmetric ansatz for the different fields is used in order 

to reduce the classical equations. The system of non-linear partial differential equations 
£NJ , obtained becomes a boundary value problem by supplementing a suitable set of boundary 

conditions. We obtain numerical evidences that the angular excitations of uncharged Q-balls, 
, which exist in flat space-time, get continuously deformed by the Maxwell and the Einstein 

terms. The electromagnetic and gravitating properties of several solutions, including the 
i- 1 spinning Q-balls, are emphasized. 

cr 

1 Introduction 

In the quest for new theoretical objects which could constitute a fraction of the exotic compo- 
nent (see e.g. [lj) of the Universe, Q-balls [21 El 0] and their gravitating counterparts, boson 
Q^ . stars El [7], emerge as serious candidates. Many efforts are carried out for an experimental 

detection of boson stars [8], but no evidence of their presence has been obtained so far. From 
the phenomenological point of view, the models containing the basic ingredients for Q-balls (i.e. 
| scalar fields) are looked for into the supersymmetric extensions of the standard model of particle 

j> , physics [9l [TOl fTTT [T2] . The super-partners of the quarks and leptons naturally lead to scalar 

•i-H . 

^ . fields. However, there is a lot of incertitude in chosing an appropriate model. 

Far away from cosmological considerations and from these elaborated lagrangians, several 
basic features of scalar fields in interaction with gravity and/or electromagnetism can be studied 
in more elementary models. Here the basic field is a complex scalar field self-interacting through 
an appropriate potential. The spectrum of Q-balls solutions turns out to be extremely rich, 
the solutions present an harmonic time dependance characterized by a 'spectraP-parameter tu s . 
Besides the spherically symmetric (fundamental) solutions, radially excitated solution exist as 
well [13]. The solutions can be done spinning both with even and odd branches [13] under the 



*E-mail: yves.brihaye@umh.ac.be 
*E-mail:thierry.caebergs@umh. ac.be 
*E-mail:terence. delsate@umh.ac.be 



1 



parity operator. Angular excitations of Q-balls in relations with the spherical harmonics are 
discussed in [Hj. Spinning solutions in super symmetric extensions of the standard model are 
constructed in [H2, [IB] . 

Once coupling the lagrangian admitting Q-balls to gravity, the solutions get continuously 
deformed by gravity but lead to an even more involved pattern [17\ 118]. Several branches occur, 
existing on finite intervals of the parameter tu s , and terminating into cusps at a series of critical 
values uj s = u s> k, k = 1, 2, ... . 

Q-balls are intimately related to a global U(l) symmetry. The gauging of this symmetry 
naturally introduces the electromagnetic field into the lagrangian. It is therefore challenging to 
understand the pattern of solutions to both gravitation and electromagnetism. The electromag- 
netic properties of spinning Q-balls are considered in [T5], Charged, spherically symmetric boson 
stars are studied in [20] for a very specific potential of the scalar field. More generic potentials 
are proposed e.g. in [22J. To our knowledge, the effects of the electromagnetic fields on the 
angular excitations of boson star, including the important case of spinning solutions have not 
yet been addressed in the litterature. The purpose of this paper is to investigate such solutions. 

In the second section of this paper, we consider the Einstein-Maxwell Lagrangian coupled to 
the Lagrangian of a scalar field admitting the basic Q-balls as generic solutions. We present the 
axially symmetric ansatz and discuss the physical quantities characterizing the solutions. 

We sketch the classical equations and the set of relevant boundary conditions in Sect. III. 
These equations were solved by using numerical methods, several results and properties of the 
numerical solutions that we obtained are presented. The solutions can be characterized by mass, 
particle number, angular momentum, electric charge and magnetic moment. In passing, we note, 
that the configurations obtained can be seen as regularized version of the Kerr-Newman black 
holes where the horizon and the essential singularity occuring at the origin are regularazied by 
the scalar field. 

2 The model 

The starting Lagrangian decribes a complex scalar field in 3 + 1 dimensions admitting Q-balls 
solutions. We then couple the scalar field to an electromagnetic field in the standard way, 
obtaining an abelian gauge theory. The resulting Klein-Gordon-Maxwell Lagrangian is then 
coupled minimally to gravity. The full action S reads: 



(1) 



where R is the Ricci scalar, G represents Newton's constant and C m denotes the matter La- 
grangian: 

C m = -^F flv F^-D^D^*-V($) , F ia , = d lk A v -d v A ik (2) 

Here <3? denotes the complex scalar fields and D„ = (<9„ — ieA^) is the covariant derivative with 
a coupling constant e. The signature of the metric is ( — h ++). The conventional form of the 
self-interacting potential is used: 

V($) = k|$| 6 -/?|$| 4 + A|$| 2 (3) 



where k, (3, A, are positive constants. It was argued in [T3] that a $ -potential is necessary for 
classical Q-ball solutions to exist. The parameter A determines the mass of the scalar field : 
(m B ) 2 = A. 

With the purpose of comparing our results with those existing in the litterature p2J HHJ E] 
in the limit e = 0, the values 

k = 1 , (3 = 2, A = 1.1 (4) 

are used in the explicit calculations. 
The energy-momentum tensor reads: 

= (D^y(D^) + (D^y(D^) + F, ia F^g a/3 + g^C rn (5) 

The conserved Noether current j^, \x = 0, 1, 2, 3, associated to the U(l)- symmetry is just 

f = $1^$*) (6) 

leading to a conserved charge Q of the system 

Q = -i I fyf^gdrdOdip (7) 



The variation of the action [T] with respect to various degrees of freedom, namely the metric, the 
scalar and vector fields leads respectively to the Einstein equations 

Gfj, u = 8nGT^ u , (8) 

(with Tn V given by ©), the Klein-Gordon equations 

dV \ 

and the Maxwell equations 

MV^gF^) = V^g + . (10) 
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3 The Q-ball equations 



3.1 Ansatz 

We will look for stationnary, axially-symmetric solutions along the direction z. For the metric, 
we use Lewis-Papapetrou coordinate frame [T?] : 

ds 2 = _ fdt 2 + L ^ dr 2 + r 2 M 2) + r 2 gin 2 + ^)2^ (n) 

where the metric functions /, L, g and w are functions of r and 9 (N.B. we use the notation 
L instead of the conventional one I to avoid confusion in the next sections). This ansatz is 
completed with the following form for the scalar and electromagnetic fields : 

<f>(t,r,6,(p) = e iuJst+ikip (f>(r,6) , A = V(r, 6)dt + A(r, 6) sin 9 dip (12) 

depending on three new functions 4>, V, A while ui s and k are constants. The periodicity of the 
scalar field = $(<p + 2ir) requires k G Z. 

Inserting the ansatz in the equations (JH]) , ([9]) , (IIOD leads to a system of seven nondinear partial 
differential equations given in appendix. These equations depend on the spectral parameter u> s , 
on the integer k and on the coupling constants e, G. It is convenient to use dimensionless 
variables and functions. In this purpose, we set 

f = M v r , Qj s = ^- , 4> = 4r , V = -TT , & = 4r , (13) 
p ' M p ' r M p ' M p ' M p ' v 7 

where M p is a mass scale associated with the phenomenology of the underlying model (see e.g. 
|21|.|16| for details). We then use a = SirGM^ and, for convenience, we drop the hat on the 
dimensionless quantities in the rest of the paper. The potential is fixed by ([!]). For e = 0, the 
Maxwell equations decouple, the remaining equations are those of the (uncharged) boson star. 
For k = 0, the equations for the fields w and A are solved trivially by w = 0, A = 0. 
Several limits of these equations have been studied over the recent years. 

• The case a = e = was the object of [13J. 

• The case e = was studied in (TTJ Q2] and in [H] . 

• Charged-rotating solutions (i.e. a = 0, e ^ 0, k = 1) are constructed in [19] 

• The case k = 0, corresponding to spherically symmetric charged boson star, was adressed 
in [20]. 
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In [13], it was argued that different angular excitations of (uncharged) Q-balls and boson-stars 
exist in correspondance with the symmetries of the spherical harmonic functions Y f k (9,ip), with 
— £ < k < i. The fundamental solutions correspond to t = k = 0. It is therefore natural to 
characterize the different families of solutions (and their corresponding charged generalisations) 
by using the integers £, k. 

3.2 Boundary conditions 

We are looking for regular, stationary solutions and presenting a symmetry (or an antisymmetry) 
under the reflexion z — > —z. The domain of integration can then be limited to r G [0, oo],9 £ 
[0, vr/2]. It turns out that the boundary conditions for the function depend on the discrete 
symmetries of the fields, i.e. on the numbers I, k . Here we give the ones corresponding to 
I = 0, 1 only (higher values can be found in [21]). The appropriate boundary conditions read: 

r = : d r f = , d r L = , g = 1 , w = , = 0, d r V = , A = 0. (14) 

for solutions with k ^ 0, while for k = solutions, we have <9 r ^| r= o = 0. The boundary conditions 
at infinity result from the requirement of asymptotic flatness and finite energy solutions: 

r = oo :/ = l, L = l, g = l, w = 0, = 0, V = 0, A = 0. (15) 

For 9 = the regularity of the solutions on the z-axis requires: 

9 = :d e f = 0, d e L = 0, g\ e=0 = 1 , d e w = , = , 8 e V = , A = 0. (16) 

for k / solutions, while for k = solutions, we have <9 e 0| r=o = 0. 
The conditions at 9 = it/ 2 are either given by 

= ~ :d e f = 0, d e L = 0, 8 e g = 0, d e w = , d e <t> = , 8 e V = , d e A = , (17) 

for even parity solutions, while for odd parity solutions the conditions for the scalar field functions 
read: 0| 9=vr / 2 = 0. 

3.3 Physical quantities 

The boson stars can be characterized by several physical parameters. The mass M and total 
angular momentum J of the solution can be computed from the appropriate Komar integral; 
it was shown in [T7] that these conserved charges can finally be read off from the asymptotic 
behaviour of the metric functions : 

M = lim r 2 8 r f , J = lim r 2 w . (18) 
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The total angular momentum J and the Noether charges Q of the two boson stars are related 
according to J = kQ as pointed out first in [23J. Boson stars with k = have thus vanishing 
angular momentum. The goal of this paper is to emphasize the electromagnetic properties of 
the solutions. The electric charge Q e and magnetic moment [i are given by the asymptotic decay 
of the electromagnetic potentials [19], i.e. 



V = Q- , A=--£- , for r^oo (19) 



and the gyromagnetic factor g can then be defined : fi = gQJ/(2M). 



4 Numerical results 

We have solved numerically the system of partial differential equations © and §§§ subject to 
the appropriate boundary conditions given in Section 3.2. This has been done using the PDE 
solver FIDISOL [25]. We have mapped the infinite interval of the r coordinate [0 : oo] to the 
finite compact interval [0 : 1] using the new coordinate z := r/(r + 1). We have typically used 
grid sizes of 150 points in r-direction and 70 points in 9 direction. The solutions presented below 
have relative errors of order 10 -3 or smaller. 



4.1 Non-rotating solutions 

This corresponds to k = 0. In this case, the system of differential equations reduces to a system 
of coupled ordinary differential equations. Four of the equations are solved trivially by g = 1, 
L = 1, w = and ^4 = 0, implying a null magnetic field. The functions f,(fi,V depend on 
the radial variable r only. In the uncharged case (i.e. with e = 0), the gravitating Q-balls 
corresponding to the potential (3), (4) exist on a finite interval of the parameter uj s , i.e. for 
oj s £ [ui,U2\- Several supplementary branches of solutions further develop in the region of 
uj s ~ iv a , forming a set of backbending branches stopping in cusp at several critical values of 
u> s . We expect this pattern of solutions to be preserved for e > 0. The evolution of the mass 
and of the conserved charge Q as functions of the electric charge Q e is reported on Fig. [T] for 
lo s = 0.85 and a = 0.1, a = 0.5. For these values, there is only one uncharged solution [T7] . 
Because charged spherically symmetric solutions were first adressed in [20] (although with a 
potential different from (3)), we put in this paper the main emphasis on several kind of angular 
excitations, especially to the important case of spinning solutions corrsponding with k = 1. 
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4.2 Spinning solutions 

Setting k = 1, the solutions are charged-spinning-gravitating Q-balls. It was shown in [T7] that 
several branches of uncharged spinning solutions exist once a > 0. Two branches of solutions 
exist on a large interval of the parameter tu s , but additionnal branches further develop with a 
smaller extension of the parameter lo s . The equations depend on three continuous parameters: 
a,e,uj s . Exploring these three independant directions represent a huge task which is not the 
scope of this paper. We therefore limited our investigation to co = 0.75 and a G [0.1,0.5] for 
which there are two branches of uncharged spinning boson stars |17| . Once increasing gradually 
the electric coupling constant e the solutions naturally get charged and develop both an electric 
and a magnetic fields. The two uncharged solutions form two branches, each labelled by the 
electric charge Q e . Examining first the case of weak gravity, it seems that the two branches can 
be continued for arbitrarily large values of the electric charge. The mass and angular momentum 
of the solutions is represented on Fig. [2] as functions of the electric charge Q e . The magnetic 
momentum \x and the corresponding constant e are supplemented in the window of the figure. 
Considering then the spinning solutions for large values of a leads to a different scenario. It 
turns out, indeed, that the two branches exist up to a maximal value of the electric charge, say 
Qe = QT^- Once this limit is attained the two branches coincide and it is very likely that 
no solution exist for Q e > Q™ ax . In the case uj s = 0.75, a = 0.5, we find Q™ ax 1.9. This is 
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Figure 2: Mass M and angular momentum J of the £ = k = 1 Q-ball as functions of the electric 
charge for a = 0.1 and lo s = 0.75 

illustrated by Fig. O we see that the electric charge links the branches of configuration emerging 
from the two (uncharged) solutions of |17j . As an attempt to interpret this phenomenon, we 
remember that, for a given value of lo s , the mass of the soliton decreases for increasing a. This 
suggests that the mass of the soliton becomes too weak for the attractive gravitational force to 
dominate the electrostatic repulsion inside the soliton. 

The contour plots of a typical charged-spinning solution in the p, z plane are presented in 
FigSlfor a = 0.5. It reveals that the functions 4>,w,A are concentrated in a torus lying in the 
equator plane, while /, V deviate only a little from spherically symmetric configurations. 

4.3 Angular excited boson stars 

Radial excitations of Q-balls have been emphasized for a long time [13], the argument demon- 
strating their existence is based on the analogy between the classical movement of a point 
particle into a potential and the stationnary radial equation of the spherically-symmetric bo- 
son star (interpreting the radial variable r as time). In [14] it was argued that Q-balls possess 
angular excitations as well and that the solutions should exist with an angular dependance in 
correspondance with the spherical harmonic functions Yg(0,tp). The gravitating counterparts 
of these excited solutions was adressed in [18], [21]. As pointed out in the previous sections, 
the cases £ = k = and £ = k = 1 respectively correspond to spherically symmetric solutions 




Figure 3: Mass M and angular momentum J of the £ = k = 1 Q-ball as functions of the electric 
charge Q e for a = 0.5 and u s = 0.75 



and spinning solutions. The case £ = 2, k = 1 corresponds to spinning solution with function cj) 
antisymmetric under the reflexion z — * —z which we discuss in the next section. 

In [J3] solutions corresponding to £ = 1, k = and £ = 2, k = 2 have been constructed. 
Here we devote some attention to the charge and gravitating version of the £ = 1, k = excited 
solution. The angular dependance in 8 goes roughly like tx cos 9. Accordingly it corresponds 
to an antisymmetric configuration of the boson field and \(j)\ presents maxima at two opposite 
points along the oz axis. The electromagnetic current is even since it is propotionnal to the 
product <p<p* . For k = 0, the equations for the fields w, A are trivially satisfied by w = A = 0; 
the solution is therefore non-spinning. The mass of this solution (for lo s = 0.75, a = 0.1) 
is presented as a function of the charge on Figj5j The contour plots of the solutions having 
Q e = 1 are shown on Fig. [6j The maxima of the energy density are localized in two regions 
centered about the points (x, y, z) = (0, 0, ±6.5). The functions /, V are even, revealing that the 
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Figure 4: Contour plots of the functions /, uu, 4>, V and A for the £ = k = 1 solution with a = 0.5 
Q e = 1 and oj s = 0.75 

two maxima of gravitating matter roughly coincide with two maxima of the density of electric 
charge. 

4.4 Parity-odd Charged spinning boson stars 

In the section above, we discussed the charged version of fundamental parity even spinning 
boson stars, it was noticed in |13| that spinning solutions with odd parity exist as well, their 
angular dependence is related to the Y% harmonic function, with Y% oc sin 6 cos 9. We succeeded 
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Figure 5: Mass of the i = 0, k = 1 solution as function of Q e (red line). Mass, angular momentum 
(black lines) and magnetic momemt (in the window) of the £ = 2, k = 1 solution corresponding 
to a = 0.1 and uj s = 0.75 

in constructing the charged-gravitating version of this solution. The mass M and angular 
momentum J are reported on Fig|5] as functions of the electric charge by the black lines. The 
corresponding magnetic momentum is reported in the window. The contour plot of such a 
solution is shown in FigEJ, we see in particular that the scalar field is maximal on a torus 
centered about the oz axis and located at z ~ 2.5. Correspondingly, the field is minimal on the 
mirror symmetric torus at z ~ —2.5. The electric and magnetic potentials are even under the 
reflexion z — > —z and present extrema in the two tori. The topology of the vector potential A 
suggests the structure of a double magnetic dipole. 

5 Conclusions 

Numerical arguments suggesting the existence of several families of charged gravitating, axially 
symmetric Q-balls have been obtained. The singularities that occur naturally in the case of 
Einstein-Maxwell solutions are avoided by the scalar field, so that the solutions are regular 
on space-time. We have argued that families of charged solitons exist for which the angular 
dependance of the scalar field possesses the symmetries of the spherical harmonics Y^(9,ip). 
These stationnary solutions can therefore be labelled by the integers £, k. 

Asymptotically, the space-time corresponding to the case i = k = solution approaches the 
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Figure 6: Contour plots of the functions /, (f) and V for the I = 1, k = solution for a = 0.5, 
uj s = 0.75 and Q e = 1 

Reisner-Nordstrom black hole. In the case £ = k = 1, the Q-balls create a Kerr-Neumann space- 
time asymptotically. The other angular excitations generates space-times which, asymptotically, 
present more involved angular dependances; it would be challenging to prove whether these 
solutions have (or not) black hole counterparts. 
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Figure 7: Contour plots of the functions /, uu, 4>, V and A for the I = 2, k = 1 solution for 
a = 0.1, uj s = 0.75 and Q e = 1 
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6 Appendix 



In this appendix we give the field equations in the axially symmetric ansatz. The notation 

— dr' 801 Liseu. 

The equations for the Maxwell fields read 
F(°< 2 ) {2fi) 2^4(°' 1 )/(°' 1 )sin(fl) 2^( 1 '°)/( 1 -°) sin(fl) LuPA^AJ^ sin 3 (g) 



+ 



+ 



^-»2 j- ^»3 f ^*2y3 

LcA4(^W(g) L^A^U^sin 3 ^) tuA^L^ sin(0) uM^L^sm^ 

/ 2 r 2 / 2 r Lr 3 Lr 

4(0,1)^(0,1) S i n (0) 2^(°' 1 ) cos(0) c/ 1 ' ) sin(0) A^u/ 1 ' ) sin(0) 



j^>3 j^>3 j^>2 

+ 2A^/(°' 1 ) cos(fl) _ ALoA/^ sin 2 (fl)cos((9) _ AuLgj cos(fl) _ Aa/°^> cos(fl) ^ 

^^^•3 ^*2,^->3 f j t"^ y*3 

2e*LV£g 2eLfus9 f^>V^> LaA^U^ sin 2 (fl) 

7 7 F 7 + 7v 

Lc^W^) sin 2 (fl) Lo^^U^ sin 2 (g) L^V^) L^°H^°) 

/ 2 r + / 2 + 2Lr 2 + 2L 

V^ ' 1 ) cot(fl) 2y( 1 -°) 2Au;cos(g)cot(g) 



A(°> 2 ) (2>0) A?Mf<M A^fW LaA^U^sin 2 ^) Lu; 2 A^ sin 2 (9) 

LwAWuW sm 2 (0) A<ML<M A^L^ A^- 1 ) cot(0) 

T 2 2L + ^ 

2Ae 2 L(f> 2 g AfW cot (5) ALu*/ ' 1 ) sin(2#) AL^ 1 ) cot(0) 2eLn<j) 2 csc(9)g 

7 + 7^ 2/ 2 r 2 2L^ + 7 

2^(0,1)^(0,1) sin (g) ^^(1,0)^(1,0) sin (g) Lo^°) sin(g) A esc 2 ((9) _ 

The scalar field equation reads 

<ft(°' 2 > + (2)0) _ A 2 e 2 L0cu 2 sin 2 (% + A 2 e 2 <\>g + 2Ae 2 Ly0c; sin(% 

^a2 ^*2y»2 j^>2 ^y^2^» 

2AeLn(j)LU 2 sin(9)g 2AeL<j>uju s wci(9)g 2Aeri(j)csc(9)g e 2 LV 2 4>g 
2eLnV (j)ojg 2eLV4>to s g Ln 2 4>uj 2 g 2Ln(j)ujuj s g L(j)uj 2 g LU'(4>)g 

J2 r J2 ^2 r 2 J2 r J2 ' 27 

n 2 0csc 2 (% ffff L^- )^ 1 ' ) g^cot(g) 2^ 1 -°) _ 
+ ^ 2L V 2 ~ ~°' 

16 



(21) 



(22) 



The equations for the metric fields read 



/(°- 2 ) (2j0) _ 6A/W ' 1 ) cot(fl) _ 3f 2 a(A(°V) 2 _ 3/ 2 a (A^ )) 2 _ aoj 2 (A™) 2 sm 2 (9) 

2Aaoj 2 A ( -°V sin(fl) cos(0) laujA^V^ 1 ) sin(fl) ao; 2 {A™) 2 sin 2 (g) 2aaM( 1 - >y( 1 -°> sin(fl) 

^3 ^a2 

4A 2 e 2 La^ 2 u; 2 sin 2 (6')fl' 8,4e 2 LFa</> 2 ^ sin(0)# 8AeLna^ 2 o; 2 sin(0)p 8AeLa0 2 ww s sin(0)g 

2AaujV(°'V cos(#) Ae 2 LV 2 a(f) 2 g 8eLnVa<j) 2 ujg 8eLVa<p 2 uj s g f(o,i) L (o,i) y(i,0) L (i,o) 
+ ^3 / ^ / + 2L^ + 2L 

/(^cot^) (Z^ 1 ^) 2 2 /( 1 -°) {f^f 4Ln 2 a<p 2 uj 2 g SLnatf 'uu s g) 4La<p 2 m 2 g 
r 2 fr 2 r f fr 2 fr f 

L(J^) 2 sin\e) + 2W^)sin 2 (fl) _ L (a/ 1,0 )) 2 sin 2 (fl) + _ a(^)) 2 _ ^ , (1>0) y 

fr f r 2 V / 



j r 2 j r y 

(23) 



3yj 2 / 2 acot 2 (fl) _ AW cos 2 (fl) _ Loo 2 sin 2 (fl) 
Lr 4 r 4 /r 2 



L(°' 2 ) (20) 2Lau; 2 (A^ 1 )) sin 2 (g) 4ALaa; 2 A^ 1 ) sin(g) cos(fl) 4LaaM^M^) sin(fl) 

r 2 ^ ' Tr 4 ^ ' 7^ 7^ 

2Laa; 2 (A( 1 ' )) 2 sin 2 (g) ALacvA^V^ sin(fl) AAfaA^cotje) 2fa (A^) 2 

J- y — J~ f ^->4 

2fa(A^) 2 4 A 2 e 2 L 2 acj) 2 a; 2 sin 2 (6)g | 4A 2 e 2 La</> 2 g | 8Ae 2 L 2 V a^uj sin(% 

^->2 ^"2^2 ^2 ^^2yt 

8AeL 2 nac/> 2 a; 2 sin(% 8AeL 2 ac/> 2 u;u; s sin(% 8AeLna</> 2 csc(% 4ALawF(°' 1 ) cos(0) 

^^2j^>2 ^^2^* ^2 ^^j^*3 

4e 2 L 2 y 2 a0 2 3 8eL 2 nFa0 2 w5( 8eL 2 Vacj) 2 uJsg AL 2 n 2 a4> 2 oj 2 g 8L 2 nac/) 2 uuu> s g 

^*2 ^y^2^» ^^2y»2 ^^2^* 

4LW^ 2 ff 4L 2 a£/(0)g 2La (V^ 1 )) 2 2La (W 1 - )) 2 4Ln 2 a0 2 csc 2 (% 
72 + 7 + 7^ + 7 + ^ 

2L(°' 1 )cot(g) (L^ ' 1 )) 2 3L( 1 -°) (L^'Q)) 2 2A 2 .Wcos 2 (fl) 2A 2 /acot 2 (fl) 
+ r 2 2Lr 2 + r 2L + /r 4 r 4 

(24) 
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g(°' 2 ) (2)0) _ 2AfgaA^ cot(fl) _ fga (A^) 2 _ fga (A^) 2 _ 3ga^ 2 (A^) 2 S m 2 (9) 
r 2 +9 Lr A Lr 4 Lr 2 fr A 

6Agauj 2 A^sm(9)cos{9) GgauA^VW sin(0) 3gau 2 (A^) 2 sm 2 (9) 

QgaujA^V^ sin(fl) QAgaujV^ cos(fl) _ 3gL (l^ ' 1 )) 2 sin 2 (fl) 3gLuuj^ sin 2 (0) 
/V /r^ 2f 2 r 2 f 2 r 

3gL (uM))\m 2 (0) g {f^ f g {f^ f 3ga (V^) 2 
2p + 2f 2 r 2 + 2p Jr 2 

3ga{V^) 2 (ff^ 1 )) 2 gd.o) (gggf 2 g L(°- 1 ) cot(g) 2gL( 1 -°) 
/ gr 2 r (7 Lr 2 Lr 

_ g(L^)) 2 _ g(L(^)) 2 2gq (^ ' 1 )) 2 / (1 0) x 2 2A 2 e 2 g 2 La^ 2 ^\9) 

2L 2 r 2 2L 2 + r 2 +zgay<f> ' j + f2r2 

6A 2 e 2 g 2 a<p 2 A 2 fgacot 2 (9) 3A 2 gauj 2 cos 2 (9) AAe 2 g 2 LVa(f> 2 uJsm(9) AAeg 2 Lna(j) 2 uj 2 sin(0) 

AAeg 2 Lacj) 2 u;u; s sm(9) YlAeg 2 na<\) 2 csc(9) 2e 2 g 2 LV 2 acj) 2 Aeg 2 LnVa(p 2 co Aeg 2 LVa<p 2 u s 

2g 2 Ln 2 oi(j) 2 uj 2 Ag 2 Lna(j) 2 i J JUJ s 2g 2 Lacf 2 oj 2 s 3gLuj 2 sm 2 (9) 2g 2 LaU(<p) 6g 2 n 2 a(f> 2 csc 2 {9) _ 
^ j2 r 2 ' ~pr ' J 2 2 f 2 r 2 J r 2 ~ 

(25) 



lo^ 2 ) m _ 8AfaivA(°V cot(g) _ Afaco (A^ 1 )) 2 AfaA^V^ csc(g) _ Afaco (A^ )) 2 
r 2 Lr A Lr A Lr 3 Lr 2 

AfaA^V^ csc(g) AA 2 e 2 a(j) 2 ojg AAe 2 Vacp 2 csc(9)g 8Aena<p 2 u; csc(fl)gg(r, 9) 

Lr r 2 r r 2 

AAeacfuJs csc(% AAfaV (0 ^ cot(g) csc(fl) AenVacj) 2 csc 2 (% 2/(°' 1 W°' 1 ) 

| 2u;/( 1 '°) 2/( 1 -°U 1 -°> An 2 atfuj esc 2 (g) g 4na0 2 a; 5 csc 2 (% | 3L^U°j 
frf r 2 r 2Lr 2 

3ujL^ 3L(i.0Ui,0) 3a; (o,i) cot(fl) 2a/ 1 - ) 4^ 2 /acj cot 2 (6Q 2a; 
2Lr + 2L + r 2 + r Lr A r 2 ~ ' 

(26) 
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